时 区 人 研究 


ARID ZONE RESEARCH 


doi:10.13866/j. azr. 2019. 02.07 


新 疆 有 无 人 为 干扰 下 土壤 盐分 估算 的 比较 


， ES. 


叶 红 云 !， 张 ” 芳 : 


(1. 新 疆 大 学 资源 与 环境 科学 学 院 ,绿洲 生态 教育 部 重点 实验 室 ,新 疆 乌鲁木齐 ”830046 ; 


2. 北京 联合 大 学 应 


文理 学 院 ,北京 


100083 ) 


摘 要 : 以 新 疆 有 无 人 为 干扰 土壤 的 实测 光谱 为 数据 源 ,利用 去 包 络 线 法 提取 其 吸收 特征 和 敏感 波段 ,建立 了 多 种 


基于 去 包 络 线 法 的 土壤 盐分 估算 的 一 元 、 多 元 线性 回归 模型 。 结 果 表 明 :Q@ 去 包 络 线 法 能 明显 增强 土壤 光谱 在 
400 ~610 nm 和 610 ~800 nm 波段 的 吸收 特征 ;@@ 有 /无 人 为 干扰 区 土壤 盐分 含量 与 去 包 络 线 一 阶 微分 的 相关 性 较 
好 ,敏感 波段 在 无 人 为 干扰 区 集中 于 可 见 光 区 (400 ~ 800 nm) ,而 人 为 干扰 区 为 可 见 光 - 近 红 外 区 (400 ~2 400 
nm); 土壤 盐分 值 预测 模型 精度 从 高 到 低 在 有 /无 人 为 干扰 区 基本 上 是 :去 包 络 线 6 个 吸收 特征 模型 > 去 包 络 线 


和 去 包 络 线 一 阶 微分 敏感 波段 模型 > 去 包 络 线 3 个 吸收 特征 模型 > 单个 吸收 特征 模型 。 有 /无 人 为 干扰 区 最 佳 预 
测 模型 的 尽 分 别 为 0.86 和 0.92。 无 人 干扰 区 最 优 模 型 (6 个 吸收 特征 ) 的 相对 分 析 误差 (RPD ) 达到 了 3. 06 ,大 于 


3 ,说 明 其 能 极 好 并 且 稳 定 的 预测 此 区 土壤 盐分 值 。 
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土壤 盐 溃 化 是 干旱 半 干 旱地 区 农业 发 展 的 主要 
障碍 因子 ,人 为 干扰 活动 又 是 导致 次 生 盐 淡化 的 直 
接 原因 。 因 此 ,准确 掌握 人 为 干扰 下 的 土壤 盐分 信 
息 ,对 于 治理 盐 渍 土 .防止 土壤 次 生 盐 涡 化 的 发 生 及 
实现 干旱 区 的 可 持续 发 展 具有 重要 意义 "”。 

光谱 分 析 技 术 是 综合 光谱 学 、 化 学 计量 学 和 计 
算 机 应 用 等 学 科 的 现代 分 析 技 术 ,与 传统 分 析 技 术 
相 比 具有 快速 无损 .无 污染 的 优点 。 当 前 ,高光 
谱 是 对 土壤 盐 渍 化 遥感 监测 的 理想 手段 之 一 ” 。 
内 外 学 者 利用 高 光谱 数据 ,以 不 同 区 域 的 土壤 盐 
分 为 对 象 展 开 了 研究 。Abrol Dehaan 和 Christian 
等 “分 别 分 析 了 土壤 高 光谱 的 特征 ,并 找到 了 盐 
分 的 光谱 敏感 波段 。 刘 迁 迁 E2929 80 O 分 别 对 
新 疆 不 同 地 区 的 土壤 盐分 进行 了 定量 分 析 。 李 上 晓 
明 . 段 鹏 程 等 "'"-… 利 用 土壤 高 光谱 盐分 敏感 波段 对 
盐分 含量 进行 了 反 演 ,发现 光谱 的 不 同 变换 形式 对 
盐分 预测 有 不 同 影响 。 彭 杰 等 …) 利用 光谱 去 包 络 
线 法 对 新 疆 南部 土壤 盐 渍 化 预测 的 效果 较 好 。 对 于 
受到 人 类 不 同 干扰 的 土壤 盐分 精准 估算 用 何 种 去 包 
络 线 法 更 为 合适 ,是 当前 需要 重点 探讨 的 问题 之 一 。 

本 研究 采用 去 包 络 线 法 提取 有 无 人 为 干扰 土壤 
的 高 光谱 敏感 波段 和 吸收 特征 ,建立 了 多 种 基于 去 
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包 络 线 法 的 土壤 盐分 估算 模型 并 检验 其 预测 效果 ， 
比较 其 对 提取 不 同 干扰 程度 下 土壤 盐 溃 化 信息 的 有 
效 性 ,提出 不 同人 为 干扰 下 土壤 盐 溃 化 的 适宜 快速 
光谱 监测 模型 。 


1 材料 与 方法 


1.1 研究 区 概况 

研究 区 在 新 疆 天 山北 荐 时 康 境内 ,地 处 山 前 冲 
R 洪 积 而 下 缘 , 地 势 较为 平坦 ,土壤 环境 背景 基本 
一 致 。 原 始 地 表土 壤 质 地 为 沙土 ,地 表 盐 渍 化 现象 
严重 ,在 西北 干旱 区 具有 代表 性 。 地 表 植 被 覆盖 度 
较 低 ,主要 为 灌木 FER., 

随 着 社会 经 济 快 速 发 展 和 人 口 的 增加 ,人 类 对 
该 区 土地 进行 了 不 同方 式 的 利用 ,对 土壤 造成 不 同 
程度 的 干扰 ,因此 将 研究 区 按 人 类 活动 干扰 程度 的 
大 小 ,分 为 无 人 干扰 区 和 人 为 干扰 区 。 无 人 干扰 区 
距 人 类 活动 区 较 远 ,地 表 基 本 保持 原始 风貌 ,植被 覆 
盖 度 较 高 ,部 分 地 区 有 成 片 梭 梭 林 ( Haloxylon am- 
modendron (C. A. Mey.) Bunge) ,伴生 植物 种 有 红 
Fil ( Tamarix juniperina Lour. ) ,SEJfUK ( Kalidium foli- 
atum (Pall. ) Moq. ) 等 。 由 于 较 少 受到 人 类 干扰 ， 
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土壤 表面 常见 黑色 生物 结 皮 ,发育 良好 ,地 衣 LIE ÉE 
结 皮 对 土壤 表面 起 着 稳定 作用 "”。 人 为 干扰 区 靠 
近 新 疆 生 产 建设 兵团 102 团 ,北部 为 弃 耕 荒地 ,植被 
TA as BE PURIA FE R EXE S 48 SERRER TE Th s PA 
部 为 人 工 林 场 ,种 有 榆树 (Ulmus pumila L. ) .枸杞 
( Lycium chinense Mill. ) 等 ,采取 滴灌 措施 。 由 于 两 
区 距离 较 近 , 仅 被 人 工 设施 分 离 , 因 此 本 研究 土壤 的 
成 士 母 质 HE 气温、 降水 等 自然 条 件 一 致 ,土壤 盐 
分 主要 受 人 类 活动 干扰 的 影响 。 
1.2 样品 采集 与 处 理 

采用 网 格 采样 法 分 别 在 无 /有 人 为 干扰 区 由 南 
向 北 布设 $ ~6 条 相距 800 ~ 1 000 m 的 采样 线 , 每 
条 线 上 选取 5 个 距离 200 ~500 m 的 能 代表 周围 环 
席 背 景 的 样 点 ,并 使 用 GPS 定位 (图 1)。 两 区 共 55 
个 样 点 ,无 人 为 干扰 区 25 个 ,人 为 干扰 区 30 个 。 
干旱 区 土壤 水 分 蒸发 强烈 ,土壤 中 的 盐分 易 随 水 分 
迁移 到 地 表 , 故 选取 能 够 反映 该 点 环境 平均 状况 的 
土壤 表层 (0 ~ 10 em) ,用 梅花 桩 采样 法 采集 每 个 样 
点 土壤 样品 1 kg ,混合 均匀 后 装 入 密封 袋 内 并 标记 
编号 。 

对 55 个 土壤 样品 经 清除 小 石 块 及 植物 根系 等 
杂 物 后 ,于 实验 室内 自然 阴干 ,随后 对 阴干 样品 过 2 
mm 第 ,将 过 筛 后 的 土壤 样品 置 于 新 样品 袋 中 ,标记 
相应 编号 送 至 中 国 科 学 院 新 疆 生 态 与 地 理 研 究 所 中 
心 试验 室 进 行 土壤 盐分 测定 。 
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图 1 样 点 示意 图 


Fig.1 Distribution of soil sampling sites 
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1.3 样品 统计 特征 

土壤 样品 含 盐 量 较 高 且 盐 分 变化 范围 较 大 ( 表 
1) ,多 为 中 度 (土壤 盐分 6~ 10 g- kg ^) .重度 盐 涡 
化 土 ( 土 壤 盐分 10 ~ 20 g- kg) db E CHE HERE 
分 >20 g- kg )5 。 从 均值 来 看 ,土壤 盐分 含量 
人 为 干扰 区 高 于 无 人 为 干扰 区 ,这 可 能 是 由 于 人 类 
活动 使 地 下 水 位 抬升 ,土壤 水 中 的 盐分 更 容易 通过 
毛细 管 上 升 至 地 表 ,使 盐分 表 聚 ,造成 了 人 为 干扰 区 
盐分 含量 较 高 "5 。 两 区 样品 数据 的 变异 系数 较 大 ， 
为 42.26% ~58.06% ,变异 性 为 中 等 变异 (10% ~ 
100% ) ,表明 样本 数据 离散 程度 较 高 ,这 有 利于 模 
型 的 构建 ,所 建 模 型 也 更 具 普 适 性 。 同 时 ,无论 是 建 
模 集 检验 集 还 是 全 部 样本 ,人 为 干扰 区 的 样品 标准 
差 和 变异 系数 均 大 于 无 人 干扰 区 。 


表 1 土壤 含 盐 量 统计 特征 
Tab.1 Statistic characteristics of soil salt 
content /(g* kg!) 
无 人 为 干扰 区 人 为 干扰 区 
建 模 集 ”检验 集 全 部 样本 建 模 集 ”检验 集 全 部 样本 
(n=15) (n=10) (n=25) (n=15) (n=10) (n=25) 
最 小 值 4.20 5.15 4.20 3.50 1.58 1.58 
最 大 值 30.08 27. 13 30. 08 41.63 36.42 41.63 
均值 16.33 15. 48 15.99 21.20 19.12 20.51 
标准 差 7.54 6.54 7.03 10.83 11.10 10.78 


数 /% 


46.17 42.26 43.95 51.09 58.06 52.54 


Bon 指 代 样品 数 。 
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1.4 光谱 采集 

土壤 光谱 数据 采用 美国 ASD 公司 的 FieldSpec 
3 型 光谱 仪 在 野外 实地 测 得 ,波长 为 350 ~ 2 500 
nm, 采样 时 间 为 2017 年 5 H 11—24 日 的 11:00 一 
16:00( 当地 时 间 ) ,天 气 晴 朗 少 云 , 基 本 无 风 。 光 谱 
测定 时 ,传感器 探头 位 于 土壤 表面 垂直 上 方 15 cm 
处 ,采用 25" 视 场 角 探头 ,每 次 采样 前 先进 行 白板 校 
正 以 去 除 暗 电流 的 影响 。 每 个 样 点 的 光谱 来 自 采 样 
点 周围 1 m 内 土壤 背景 相近 的 5 处 土壤 ,每 处 采集 
1 次 ,每 次 10 条 光谱 , 共 50 条 。 
1.5 光谱 预 处 理 及 吸收 特征 参数 提取 

采用 ASD View Spec Pro 对 光谱 预 处 理 , 运 用 
Origin Pro 9. 1 对 光谱 反射 率 进 行 Savitaky-Golay 平 
滑 去 噪 处 理 。 由 于 土壤 反射 光谱 前 段 350 ~ 399 nm 
和 后 段 2 401 ~2 500 nm 的 噪声 较 大 ,因此 剔除 这 两 
段 光谱 。 此 外 ,受到 水 分 吸收 带 的 影响 ,在 1 400 
nm,1 900 nm 和 2 200 nm 处 的 波动 也 较 大 " ,因此 


剔除 反射 率 值 波动 剧烈 的 1359 ~1 410 nm,1 820 ~ 
1 942 nm 和 2 197 ~2 210 nm 波段 。 利 用 Envi 5.1 
进行 去 包 络 线 人 处 理 后 再 用 Origin Pro 9. 1 对 其 进行 
微分 处 理 。 

去 包 络 线 变 换 后 ,提取 如 下 几 个 典型 吸收 特征 : 
最 大 吸收 深度 、 吸 收 峰 总 面积 和 最 大 吸收 宽度 ( 
2 ) 用 以 建 模 "…” 。 
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图 2 光谱 吸收 峰 吸收 特征 示意 图 


Fig.2 Schematic diagram of spectral absorption features 


最 大 吸收 深度 (D) 即 吸收 峰 内 去 包 络 线 最 低 点 到 

归 一 化 包 络 线 间 的 垂直 距离 。 公 式 为 : 
Dzl-p, 
式 中 :pu 为 吸收 峰 内 去 包 络 线 最 低 点 光谱 去 包 络 线 
值 。 
吸收 峰 总 面积 (4) 即 去 包 络 线 值 积分 ,是 指 吸 

收 波长 范围 内 光谱 曲线 与 去 包 络 线 所 围 成 的 面积 
小 。 公 式 为 : 


4 = 三 1 -f(w)dw 


式 中 :f(w) 即 为 光谱 曲线 ;wi wa 为 吸收 峰 两 肩 波长 
值 。 

最 大 吸收 宽度 ( WW) 即 为 吸收 峰 两 肩 之 间 的 距 

离 。 公 式 为 : 
Wzw,-w, 
模型 构建 与 精度 评价 
采用 Kennard-Stone 法 将 两 种 人 为 干扰 的 样品 
盐分 数据 划分 为 建 模 集 和 检验 集 ,无 人 为 干扰 区 的 
建 模 集 样本 数 为 15 ,检验 集 为 10, 人 为 干扰 区 的 个 
数 分 别 是 20 和 10( 表 1)。 

建 模 方法 采用 一 元 线性 回归 和 多 元 线性 回归 ， 
模型 精度 采用 决定 系数 (RR )、 均 方 根 误差 (root 
mean square error, RMSE ) ,相对 分 析 误 差 (relative 
percent deviation, RPD) 3 种 指标 来 评价 ,其 计算 公 
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式 如 下 。 
2 Èra (x) 
R= p 

22 sa Cg -x) 


I n 2 
RMSE aJ ius =k) 


JEE C 2) 

RED = RMSE 
XB in 表示 样本 个 数 ;x, 表示 模型 预测 值 ;xo 表示 
实测 值 ;x 表示 样品 平均 值 。 

R 反映 预测 值 解释 实测 值 变 化 的 程度 , RMSE 
用 来 检验 模型 模拟 值 和 实测 值 的 符合 度 。RPD 用 
来 验证 模型 的 稳定 性 和 预测 能 力 。 若 R 越 大 ， 
RMSE 越 小 ,说 明 模型 精度 越 高 。 此 外 , 当 1. 5 < 
RPD «2. 0 时 ,表明 模型 只 能 对 样品 高 含量 与 低 含 
量 进 行 粗略 估计 ,2. 0 RPD «2.5 时 ,表明 模型 具 
有 较 好 的 定量 预测 能 力 ;2.5 RPD «3. 0 时 ,模型 
具有 很 好 的 预测 能 力 ;RPD=3. 0 时 ,分 析 模 型 具有 
极 好 的 预测 能 力 "*)。 


2 结果 分 析 


2.1 有 无 人 为 干扰 下 土壤 高 光谱 特征 

分 别 平均 无 人 为 干扰 区 、 人 为 干扰 区 的 光谱 后 ， 
得 到 2 条 原始 反射 率 曲线 (图 3a) 。 从 图 中 可 以 看 
出 ,两 区 反射 率 曲线 走势 基本 一 致 , 在 可 见 光 范围 内 
580 nm 处 有 交叉 节点 ,节点 前 无 人 为 干扰 区 土壤 反 
射 率 高 于 人 为 干扰 区 ,因为 无 人 为 干扰 区 盐分 均值 
低 于 人 为 干扰 区 ,表现 出 盐分 低 反 射 率 高 的 规律 。 
而 节点 之 后 ,人 为 干扰 区 的 反射 率 较 无 人 为 干扰 区 
的 高 ,呈现 出 盐分 高 反射 率 高 的 特征 , 且 随 着 波长 的 
增加 ,两 区 反射 率 值 差异 随 之 增 大 。 男 外 ,曲线 在 
1 416 nm,1 947 nm 和 2 216 nm 波段 处 有 明显 的 吸 
收 峰 。 

两 区 光谱 的 去 包 络 线 (CR) 曲线 (图 3a) 中 
1 416 nm 1 947 nm 和 2 216 nm 波段 靠近 水 分 吸收 
带 处 的 吸收 特征 得 到 加 强 , 变 得 尖锐 且 明 显 。 而 可 
见 光 范 围 内 ,两 区 在 400 ~ 610 nm 和 610 ~ 800 nm 
波段 处 从 原始 光谱 反射 率 曲线 中 很 难 发 现 的 微弱 吸 
收 特征 得 到 放大 , 且 CR 值 差异 较 大 。 这 表明 CR 法 
在 可 见 光 范围 内 对 土壤 光谱 的 反射 和 吸收 特征 放大 
效果 较 好 ,尤其 是 突出 了 不 同人 为 干扰 下 土壤 光谱 
的 吸收 差异 ,因此 可 利用 此 处 的 差异 分 析 有 无 人 为 
干扰 的 土壤 盐分 值 。 
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注 :CR 表示 去 包 络 线 ,FDCR 表示 去 包 络 线 一 阶 微分 。 下 同 。 


图 3 ”有 无 干扰 下 土壤 原始 光谱 反射 率 去 包 络 线 曲线 (a) 和 去 包 络 线 一 阶 微分 (b) 


Fig.3 Spectral curves and continuum removal (a) and continuum removal first — order differential (b) of soils under 


different intensities of human disturbance 


两 区 光谱 去 包 络 线 经 一 阶 微分 (FDCR ) 变换 之 
后 的 曲线 (图 3b) 在 大 部 分 波段 都 较 接近 , 即 FDCR 
值 差异 不 大 。400 ~ 800 nm 内 有 几 处 较 明 显 的 差 
异 , 分 别 为 430 nm,572 nm .664 nm 和 692 nm, 而 到 
近 红 外 波段 范围 内 , 仅 有 1 356 nm 和 2 400 nm 附近 
有 所 区 分 。 

将 CR 差异 明显 的 2 个 吸收 峰 放大 得 到 ,400 ~ 
610 nm 波段 的 吸收 峰 较 平滑 (图 4a) ,CR 值 在 无 人 
为 干扰 区 的 最 小 值 为 0.937(498 nm) ,人 为 干扰 区 
为 0.912(505 nm) 。 此 波段 无 人 为 干扰 区 的 CR 值 
大 于 人 为 干扰 区 , 即 吸 收 深度 和 吸收 面积 均 是 无 人 
为 干扰 区 较 小 ,又 由 于 人 为 干扰 区 土壤 盐分 值 高 于 
无 人 干扰 区 ,因此 具有 随 着 盐分 的 增加 而 增 大 的 吸 
收 特征 。 

而 到 红 光 附近 (610 ~ 800 nm) 吸收 峰 则 较 尖 锐 
突出 (图 4b) ,尤其 无 人 为 干扰 区 的 CR 值 在 681 nm 
波段 处 达到 了 0.985 , 比 人 为 干扰 区 在 683 nm 处 的 
最 低 值 0.993 低 0.012。 与 前 一 个 吸收 峰 的 吸收 规 
律 相反 ,此 波段 的 光谱 CR 值 是 无 人 为 干扰 区 高 于 
人 为 干扰 区 ,呈现 出 盐分 高 吸收 弱 的 特点 。 
一 一 无 人 为 干扰 区 
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而 到 红 光 附近 (610 ~ 800 nm) 吸 收 峰 则 较 尖 锐 
突出 (图 4b) ,尤其 无 人 为 干扰 区 的 CR 值 在 681 nm 
波段 处 达到 了 0.985 , 比 人 为 干扰 区 在 683 nm 处 的 
最 低 值 0.993 低 0.012。 与 前 一 个 吸收 峰 的 吸收 规 
律 相 反 ,此 波段 的 光谱 CR 值 是 无 人 为 干扰 区 高 于 
人 为 干扰 区 ,呈现 出 盐分 高 吸收 弱 的 特点 。 

2.2 有 无 人 为 干扰 下 土壤 含 盐 量 与 高 光谱 的 相关 
性 分 析 

无 人 为 干扰 区 土壤 含 盐 量 与 CR 在 599 ~ 711 
nm,1 358 nm,1 817 ~1 992 nm,2 288 ~2 313 nm, 
2 354 nm 和 2 387 nm 波段 处 为 正 相关 ,其 余 为 负 相 
关 ( 图 Sa) 。 最 大 相关 系数 分 别 在 677 nm 2 229 nm 
波段 ,为 0.62 和 0.69。 此 区 土壤 盐分 含量 与 FDCR 
的 相关 性 较为 复杂 , 正 负 相 间 ,缺乏 连续 性 ,但 由 于 
微分 处 理 能 够 将 曲线 中 一 些 隐 上 临 的 信息 放大 ,使 得 
许多 波段 的 相关 性 明显 增强 。657 nm,664 nm 和 
693 nm 处 的 相关 系数 值 排 在 前 3 位 ,最 大 值 达 到 

-0.70, 因 此 反映 出 无 人 为 干扰 区 土壤 光谱 FDCR 
的 敏感 波段 基本 处 在 可 见 光 范 围 内 。 
人 为 干扰 区 的 CR 与 盐分 含量 的 相关 性 波动 较 


去 包 络 线 值 
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图 4 有 无 人 为 干扰 区 土壤 光谱 在 不 同 波段 的 去 包 络 线 值 


Fig.4 Values of continuum removal in different bands under different intensities of human disturbance 
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有 无 人 为 干扰 下 土壤 含 盐 量 与 光谱 去 包 络 线 及 其 一 阶 微分 的 相关 性 


Fig.5 Correlation between soil salt content and continuum removal and its first — order differential under different intensities 


of human disturbance 


大 (图 5b)。 正 相关 的 波段 变 多 , 且 相 关系 数 也 更 

可 见 光 、 近 红外 两 部 分 内 的 相关 系数 最 大 值 分 
别 在 678 nm 和 1 944 nm 波段 取得 ,为 0. 49 和 
0.56 ,与 无 人 为 干扰 区 类 似 ,CR 通过 一 阶 微分 处 理 
后 与 盐分 的 相关 系数 有 明显 提高 ,人 为 干扰 使 得 土 
二 光谱 FDCR 的 敏感 波段 在 可 见 光 - 近 红 外 波段 内 
均 有 分 布 。 最 大 相关 系数 在 677 nm、2 010 nm 和 


1 549 nm 处 ,分别 为 0.62、- 0.61 和 -0.58, 比 CR 
与 盐分 的 最 大 相关 系数 所 高 了 6% 。 


2.3 ”有 无 人 为 干扰 土壤 光谱 CR FDCR 的 敏感 
波段 选取 

人 为 干扰 使 土壤 盐分 含量 产生 变化 ,因此 有 无 
人 为 干扰 土壤 光谱 特征 出 现 差 异 ,导致 有 无 人 为 干 
扰 土壤 光谱 敏感 波段 有 所 不 同 。 根 据 显 著 (P = 
0.01 ) 且 相关 系数 大 的 原则 ,分别 选取 不 同人 为 干 
扰 土壤 光谱 CR FDCR 的 敏感 波段 。 由 于 在 下 述 提 
取 CR 吸收 特征 时 是 在 2 个 波段 共 提 取 6 个 吸收 特 
征 ,为 了 方便 对 比 模型 精度 ,敏感 波段 也 选择 6 个 。 

在 无 人 为 干扰 区 CR 与 盐分 相关 系数 通过 0.01 
水 平 显著 性 检验 的 2 个 波段 处 各 选择 3 个 敏感 波 
段 ,644 ~ 694 nm 波段 内 选 653 nm,668 nm 和 677 
nm,2 190 ~2 246 nm 内 为 2 229 nm、2 230 nm 和 
2 244 nm, HTF FDCR 的 敏感 波段 多 分 布 在 可 见 光 
范围 内 ,因此 在 可 见 光 范 围 内 选取 657 nm ,664. nm, 
675 nm .692 nm,693 nm 和 699 nm, 

人 为 干扰 区 CR 的 敏感 波段 选取 相关 系数 排 
在 前 列 的 678 nm、679 nm,682 nm, 686 nm, 
1 944 nm 和 1 949 nm 波段 。FDCR 的 敏感 波段 则 
分 别 为 477 nm .661 nm,676 nm .677 nm 1 549 nm 
和 2 010 nm, 


2.4 土壤 含 盐 量 高 光谱 定量 反 演 模型 的 构建 与 验证 
2.4.1 基于 CR 及 FDCR 敏感 波段 的 估算 模型 
为 了 考查 CR 及 FDCR 敏感 波段 建 模 对 有 无 人 为 干 
扰 土壤 盐分 的 预测 效果 ,分 别 选取 CR 和 FDCR 
敏感 波段 对 盐分 含量 构建 多 元 线性 模型 ( 表 2) 。 

CR 敏感 波段 为 自 变量 的 模型 中 ， npud 
型 精度 优 于 人 为 干扰 区 ,其 RO 和 RMSE 分 别 为 
0.57 和 4.79 ,说 明 CR 法 对 无 人 为 干扰 区 的 土壤 光 
谱 信 息 放 大 效果 较 好 。 


表 2 有 无 人 为 干扰 下 光谱 去 包 络 线 及 其 一 阶 微 

分 敏感 波段 模型 
Tab.2 Inversion models for soil salt content on sensitive 
bands of continuum removal and its first-order differential 


under different intensities of human disturbance 


无 人 干扰 区 人 为 干扰 区 
R? RMSE R? RMSE 
CR 0.57 4.79 0.46 7.74 
FDCR 0.71 3.94 0.67 4.10 


注 :R 表示 决定 系数 ,RMSE 表示 均 方 根 误差 。 下 同 。 


光谱 去 包 络 线 经 一 阶 微分 变换 后 ,利用 其 敏感 
波段 所 建 模型 精度 有 较 大 的 提高 ,尤其 人 为 干扰 区 ， 
R? 提高 了 0.21,RMSE 降低 了 3.64。 这 说 明 微 分 处 
A ee a 能 够 放大 原 
本 的 曲线 信息 ,提高 一 些 波段 的 相关 性 。 虽 然 人 为 
干扰 区 土 拉 受到 人 类 沥 和 、 FEK, PR 
理 和 化 学 性 质 都 有 所 变化 ,使 得 土壤 光谱 也 变 得 更 
T. 但 一 阶 微分 能 ee i 从 

提高 光谱 与 盐分 的 相关 性 。 因 此 去 包 络 线 一 阶 微 
2 uu nt re 模型 精 
度 具 有 关键 作用 。 


T 旱 区 
2.4.2 基于 CR 吸收 特征 的 估算 模型 ”计算 出 每 


个 土 样 光谱 去 包 络 线 在 400 ~ 610 nm 和 610 ~ 800 
nm 两 个 波段 的 最 大 吸收 深度 .吸收 面积 和 最 大 吸收 
宽度 。 先 以 每 种 吸收 特征 参数 作为 单一 自 变 量 构建 
盐分 拟 合 的 一 元 一 次 线性 模型 ,然后 把 上 述 2 个 波 
段 内 每 个 波段 的 3 个 吸收 特征 综合 起 来 进行 建 模 ， 
最 后 将 2 个 波段 内 共 6 种 吸收 特征 参数 一 起 作为 自 
变量 建立 多 元 线性 模型 来 预测 盐分 含量 ( 表 3)。 

单 因 素 模 型 中 ,无 论 是 在 哪个 区 哪个 波段 ,吸收 
面积 和 吸收 深度 的 RR 均 大 于 吸收 宽度 的 R , RMSE 
则 相反 。 其 中 ,以 无 人 为 干扰 区 400 ~ 610 nm 波段 
的 吸收 面积 为 自 变 量 的 模型 精度 最 高 , R 达到 
0.72, 其 次 是 此 区 同 波段 的 吸收 深度 模型 ,R 为 
0.70 ,随后 是 人 为 干扰 区 400 ~ 610 nm 波段 的 吸收 
深度 .吸收 面积 模型 R 分别 为 0.65 ,0.64。 总 体 来 
看 , 除 人 为 干扰 区 的 吸收 深度 模型 外 ,不 同 程 度 人 为 


干扰 下 的 土壤 盐分 反 演 模型 均 是 400 ~610 nm 波段 
的 优 于 610 ~ 800 nm 波段 。 

多 因素 建 模 中 ,对 比 2 波段 的 综合 吸收 特征 模 
型 发 现 ,有 无 人 为 干扰 区 均 是 400 ~ 610 nm 波段 的 
模型 精度 高 于 610 ~ 800 nm 波段 。 但 400 ~ 610 nm 
波段 的 模型 是 无 人 为 干扰 区 优 于 人 为 干扰 区 , 而 
610 ~800 nm 波段 的 则 相反 ,说 明 不 同人 为 干扰 土 
壤 光 谱 在 两 波段 内 的 吸收 特征 响应 有 所 不 同 。2 波 
段 6 个 吸收 特征 综合 建 模 无 人 干扰 区 精度 较 高 ,R 
达到 0.85 , 相 比 单 因素 建 模 有 显著 提高 。 说 明 吸 收 
深度 .面积 和 宽度 仅 能 部 分 体现 光谱 的 吸收 特征 ,不 
够 全 面 ,而 综合 2 波段 共 6 种 吸收 特征 则 能 较 好 地 
反映 出 土壤 的 光谱 吸收 状况 。 
2.4.3 模型 的 验证 与 最 佳 预测 模型 选择 ”利用 独立 
于 建 模样 本 的 检验 样本 对 敏感 波段 模型 和 吸收 特征 
模型 进行 对 比 验证 ,以 挑选 出 最 佳 预测 模型 ( 表 4) 。 


RI ”有 无 人 为 干扰 下 盐分 含量 吸收 特征 参数 模型 


Tab.3 [Inversion models for soil salt content based on spectral absorption index under different intensities 


of human disturbance 


9 变量 无 人 为 干扰 区 人 为 干扰 区 
光谱 吸收 参数 R RMSE R? RMSE 
单 因素 400 ~610 nm 吸收 深度 0.70 3.98 0.65 6.19 
吸收 面积 0.72 3.63 0. 64 6.31 
吸收 宽度 0. 46 5.34 0.30 8.78 
610 ~ 800 nm 吸收 深度 0.44 5.43 0.46 7.74 
吸收 面积 0.52 5.05 0.45 TT 
吸收 宽度 0.33 5.95 0.34 8.53 
多 因素 400 ~610 nm 3 个 吸收 特征 0.77 3.51 0.68 5.98 
610 ~ 800 nm 3 个 吸收 特征 0. 50 5.01 0.53 4.97 
2 波段 6 个 吸收 特征 0. 85 3.01 0. 79 4.13 
表 4 不 同 模型 预测 土壤 盐分 效果 
Tab.4 Soil salt contents predicted by different models 
m 自 变量 无 人 为 干扰 区 l 人 为 干扰 区 
R RMSE RPD R RMSE RPD 
敏感 波段 CR 0.81 2.56 2.06 0. 80 4.23 2.11 
FDCR 0.76 3.03 1.94 0.81 2.12 2.55 
吸收 特征 400 ~610 nm 吸收 深度 0.62 3.81 1.88 0.40 8.13 0.82 
吸收 面积 0.53 4.24 1.07 0.27 9.00 0. 78 
吸收 宽度 0.20 5.55 0.50 0.67 6.06 1.72 
610 ~ 800 nm 吸收 深度 0.46 4.54 0.93 0. 68 4.98 1.85 
吸收 面积 0.40 4.82 0.81 0.33 8.65 0.70 
吸收 宽度 0.27 5.30 0.61 0.18 10.10 0.30 
400 ~610 nm 3 个 吸收 特征 0.78 4.58 2.76 0.62 3.79 1.30 
610 ~ 800 nm 3 个 吸收 特征 0.47 4.53 0.93 0. 60 5.78 1.52 
2 波段 6 个 吸收 特征 0.92 1.78 3.06 0. 86 2.11 2.93 


注 :RPD 表示 相对 分 析 误 差 。 


在 所 有 模型 中 ,6 种 吸收 特征 共同 作为 自 变量 建立 
的 模型 精度 高 于 其 余 模型 ,有 、 无 人 为 干扰 区 的 RR 
分 别 达 到 了 0. 86 和 0. 92。 其 中 ,后 者 模型 精度 最 
高 RMSE 值 仅 为 1.78, 是 所 有 模型 中 最 小 的 , RPD 
值 为 3. 06 ,达到 所 有 模型 中 的 最 高 值 , 且 大 于 3 ,这 
表明 利用 吸收 特征 联合 对 无 人 干扰 土壤 盐分 进行 预 
测 具 有 稳定 且 极 好 的 效果 。 人 为 干扰 区 以 6 种 吸收 
特征 为 自 变 量 建立 的 多 元 线性 模型 RPD 值 2. 93 
(>2.5) , 即 能 够 很 好 地 预测 人 为 干扰 下 的 土壤 含 
盐 量 。 

除了 以 上 两 个 模型 ,利用 CR FDCR 建立 的 预 
测 模型 精度 也 相对 较 高 ,其 R 在 有 /无 人 为 干扰 区 
分 别 达到 0. 81,0. 76 和 0.80 .0.81,RPD 分 别 达 到 
2.06 .1.94 和 2.11.2.05 ,表明 光谱 变换 的 敏感 波段 
建 模 虽然 不 是 最 佳 模型 ,但 其 稳定 性 好 ,也 能 够 较 好 
的 对 有 无 人 为 干扰 的 土壤 盐分 进行 预测 ,尤其 人 为 
干扰 区 的 敏感 波段 模型 RPD 值 都 大 于 2。2 波段 的 
3 个 吸收 特征 分 别 综合 的 模型 仅 无 人 干扰 区 400 ~ 
610 nm 波段 的 预测 效果 较 好 ,RPD >2。 而 利用 一 
种 吸收 特征 建立 的 模型 中 ,只 有 无 人 为 干扰 区 的 
400 ~ 610 nm 波段 的 吸收 深度 模型 (RR 和 RPD 分 别 
为 0.62 和 1.88) ,以 及 人 为 干扰 区 的 400 ~ 610 nm 
波段 的 吸收 宽度 模型 (R 为 0.67,RPD 为 1.72) 和 
610 ~ 800 nm 波段 的 吸收 深度 模型 (RR 达到 了 
0.68,RPD 为 1.85) 的 RPD >1.5, 能 粗略 估计 盐分 
含量 ,其 余 模型 则 较 差 。 

虽然 利用 基于 CR 的 敏感 波段 和 吸收 特征 建立 
的 多 元 线性 模型 都 能 较 好 甚至 极 好 地 预测 有 无 人 为 
干扰 土壤 的 含 盐 量 ,但 6 种 吸收 特征 联合 建立 的 模 
型 在 两 区 内 对 模型 精度 提高 效果 最 显著 ,因此 ,无 论 
有 无 人 为 干扰 ,最 佳 模 型 都 是 利用 CR 提取 出 的 吸 
收 特征 建立 的 多 元 线性 模型 。 


3 结论 

本 文 利用 CR 提取 出 的 敏感 波段 和 吸收 特征 建 
立 了 有 无 人 为 干扰 土壤 盐分 含量 的 高 光谱 估算 模 
型 ,重点 对 比 了 各 种 模型 在 无 人 干扰 和 人 为 干扰 区 
土壤 盐分 预测 精度 的 差异 。 具 体 结论 如 下 : 

(1) 光谱 曲线 进行 CR 处 理 后 ,吸收 带 显著 增 
强 ,尤其 可 见 光 范围 内 400 ~ 610 nm 和 610 ~ 800 
nm 波段 处 的 吸收 峰 被 明显 放大 。 而 对 FDCR 变换 
后 ,无 人 干扰 区 的 敏感 波段 主要 分 布 于 可 见 光 
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(400 ~ 800 nm) 范围 内 ,人 为 干扰 区 则 是 可 见 光 - 
近 红 外 波段 (800 ~2 400 nm) 。 

(2) 利用 敏感 波段 建立 的 模型 中 , 相 较 于 无 人 
干扰 区 ,人 为 干扰 区 的 FDCR 模型 精度 更 高 , RPD 
大 于 2.5, 具 有 很 好 的 预测 土壤 含 盐 量 的 能 力 。 说 
明 在 CR 的 基础 上 再 进行 一 阶 微分 处 理 能 提高 人 为 
干扰 区 土壤 盐分 预测 精度 。 因 此 ,对 于 土壤 盐分 变 
异 系数 大 ,人 类 活动 干扰 强烈 地 区 的 土壤 含 盐 量 预 
测 更 适合 用 FDCR 模型 。 

(3) 总 体 上 无 论 有 无 人 为 干扰 ,模型 精度 是 去 
包 络 线 6 个 吸收 特征 > 去 包 络 线 敏 感 波段 > 去 包 络 
线 3 个 吸收 特征 > 单个 吸收 特征 模型 。 因 此 ,基于 
去 包 络 线 的 6 种 吸收 特征 建立 的 多 元 线性 模型 都 是 
预测 精度 最 高 且 最 稳定 的 模型 , 反 演 效果 十 分 理想 ， 
无 /有 人 为 干扰 区 R 分 别 为 0.92 和 0. 86, RPD 为 
3.06 和 2.93 ,可 极 好 、 很 好 地 预测 土壤 盐分 。 而 基 
于 CR 敏感 波段 的 模型 除 无 人 为 干扰 区 的 FDCR 
外 ,其 余 模型 R 均 大 于 0.8, 且 RPD 大 于 2, 也 对 土 
壤 盐 分 具有 较 好 的 预测 能 
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Estimation and Comparison of Soil Salinity under Different Intensities 


of Human Disturbance in Xinjiang 


WANG Ning, XIONG Hei-gang, MA Li-fang , YE Hong-yun', ZHANG Fang! 
(1. College of Resources and. Environment Sciences, Xinjiang University, Urumqi 830046, Xinjiang , China ; 
2. College of Arts and Science, Beijing Union University, Beijing 100083 , China) 


Abstract: Soil spectra under different intensities of human disturbance were used as the data sources to extract 
the absorption characteristics and sensitive bands by the continuum removal method. Linear regressions were estab- 
lished for estimating soil salinity based on the continuum removal method. The results showed that; (1) The continu- 
um removal method could be used to significantly enhance the absorption characteristics of soil spectra in wave- 
length bands of 400 -610 nm and 610 -800 nm; (2) The correlation between soil salt content in the manned/un- 
manned disturbance area and the first derivative of the continuum removal ( FDCR) was better than continuum re- 
moval ( CR) , the sensitive bands were concentrated in the Vis area in the unattended area (400 — 800 nm), and 
the artificially disturbed area was Vis-NIR (400 -2 400 nm) ; (3) Basically, the accuracy of the models for predic- 
ting soil salinity in the manned/unmanned disturbance area was in an order of 6 absorption characteristic models of 
continuum removal > sensitive band model of CR and FDCR >3 absorption characteristic models of continuum re- 
moval » single absorption feature model. The values of R^ of the best predictive model of the manned/unmanned 
disturbance regions were 0. 86 and 0.92, respectively. The value of the optimal model of relative percent deviation 
for the unmanned disturbed zone reached 3. 06 and higher than 3, which revealed that the model was excellent and 
stable in predicting soil salinity in this area. 


Key words: human disturbance; spectral continuum removal; soil salinity; hyperspectrum; Xinjiang 


